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Q Abstract: A numerical study of multi-phase granular materials based upon micro- 

^ mechanical modelling is proposed. Discrete element simulations are used to inves- 

■J^ tigate capillary induced effects on the friction properties of a granular assembly in 

r^ the pendular regime. Capillary forces are described at the local scale through the 

. Young-Laplace equation and are superimposed to the standard dry particle inter- 

C^ action usually well simulated through an elastic-plastic relationship. Both effects 

Q of the pressure difference between liquid and gas phases and of the surface ten- 

sion at the interface are integrated into the interaction model. Hydraulic hysteresis 
^ is accounted for based on the possible mechanism of formation and breakage of 

P capillary menisci at contacts. In order to upscale the interparticular model, triax- 

ial loading paths are simulated on a granular assembly and the results interpreted 
through the Mohr-Coulomb criterion. The micro-mechanical approach is validated 
^ with a capillary cohesion induced at the macroscopic scale. It is shown that inter- 

■^ particular menisci contribute to the soil resistance by increasing normal forces at 

contacts. In addition, more than the capillary pressure level or the degree of satu- 
ration, our findings highlight the importance of the density number of liquid bonds 
on the overall behaviour of the material. 



CO 

O 

CN Keywords: Discrete Element Method, micromechanics, capillarity, multi-phase 
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•i-H 

K^ 1 Introduction 

C^ An outstanding feature of granular materials lies in the existence of different scales 

of interest. First, the microscopic scale corresponding to the scale of individual 
particles can be considered. The microscopic scale is intimately related to the in- 
teraction processes that occur between adjoining particles. Second, a mesoscopic 
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scale can be identified, associated to a set of a few particles. The mesoscopic scale 
is relevant to describe local kinematics ( jCambou, Chaze, and Dedecker|pOOO| )), or 



to address the creation and deletion of force chains amongst several elements (Rad 



jai', Roux, and Moreau| ( |1999) ). Third, the macroscopic scale which is the scale of 
the specimen, generally denoted as the material scale. 

As a matter of fact, the overall mechanical behaviour of granular materials is in- 
timately related to the local properties that take place at the microscopic scale. In 
addition to dry contact interactions, the presence of a surrounding liquid can mod- 
ify local deformation, attrition, and sliding between particles. The material is thus 
a tri-phasic medium where the skeleton, liquid and air interact depending on the 
thermodynamic equilibrium. According to the liquid content and to the void ra- 
tio of the medium, different saturation regimes can be observed. For small liquid 
contents, as long as the liquid tends to be concentrated as independent menisci 
between adjoining particles, this is the pendular regime. For increasing liquid con- 
tent, menisci begin to merge. From this point, all liquid-gas configurations until 
the gaseous phase consists in gas bubbles within the liquid are classified as funicu- 
lar configurations. Beyond, the medium is in a capillary state until it tends toward 
a fully saturated regime. For each case, the presence of both liquid and gas in- 
duces attractive forces that confer specific constitutive properties to the material. 
However, even though interparticular menisci have been dealt with a great atten- 



tion in the pendular regime (Hotta, Takeda, and lonya (1974); Lian, Thornton, and 



Adams| ( |T993| ); |Soulie, Cherblanc, Youssoufi, and Saix| ( [2006| )), their modelling for 



higher saturation regimes is still a great challenge to deal with due to complex ther- 



modynamic arrangements between the three phases (Urso, Lawrence, and Adams 



( |1999| ); |Murase, Mochida, and Sugama| ( |2004[ )). The present paper focuses on liq- 
uid bridges as defined in the pendular state as an attempt toward the understanding 
of induced complex macroscopic behaviours. The analysis is therefore limited to 
low saturation degrees where the pendular assumption can be assumed. 

Assessing the constitutive behaviour of such multi-phasic materials at the macro- 
scopic scale can be addressed through sophisticated phenomenological constitutive 



models, requiring specific assumptions and mathematical refinements (see Cam- 



bou (1998) for a general review). This line of thinking has prevailed during the 



past decades. A powerful alternative concerns the micromechanical approaches, 
in which the description of physical phenomena is dealt with at the microscopic 
scale. The overall behaviour is therefore derived thanks to homogeneous schemes, 
relating both microscopic and macroscopic scales as presented in |Chr istoffersen,| 
[Mehrabadi, and Nemat-Nasser] ( |1981| ); |Cambou| ( |1998| ); [Kruyt and~Rothenburgl 
( [2000| ); [Nlcot and Darve] ( [2005] ) or |Lu and Likos] ( [2006| ) for example. Such meth- 
ods are appealing since the local physics can be described by simple equations. 
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requiring a small number of parameters with generally a clear physical meaning. 

Following a micromechanical approach, discrete element methods (DEM) consti- 
tute a particularly relevant computational tool to track the constitutive response of 
granular masses along a variety of loading paths. Indeed, the rheological behaviour 



is obtained without any global hypotheses. Since the pioneering work of Cundall 



and Strack] ( |1979| ), DEM approach has been popularized in the field of granular 
materials, and extension to multi-physical coupling issues is an emergent line of 
research. For example, hydro-mechanical modelling supports now a great source 



of efforts, both for unsaturated ( [Richefeu, Youssoufi, Peyrous, and Radjai| ( |2007| ); 
Shamy and Groger| ( |2008| ); [Scholtes, Chareyre, Nicot, and Darve (2009])) and satu- 
rated flows ( |Han, Feng, and Owen| ( |2007| ); |Zeghal and Shamyl ( [2Q08) ). This paper is 
situated in this context, by showing how the local parameters characterizing wetted 
contact interactions within a granular medium can influence the overall response of 
the material. In the following, the numerical specimen is assumed to be homoge- 
neous and constituted by a sufficient number of grains so as a constitutive relation 
can be derived from both homogeneous strain and stress tensors: the specimen is 
considered as a Representative Volume Element (RVE) of a granular medium 



2 Micromechanical model 

DEM has been extensively used to study soil mechanics providing, for instance, 
some insights into shear strength and deformation properties of granular soils (see 
Iwashita and Oda| ( [l988| ) and [Ting, Corkum, Kauffman, and Greco] ( |1989| ) for ex- 
ample). Recently, DEM has been enhanced to investigate unsaturated granulates 



features by considering the possible effects of capillary water between grains ( |Gili 
and Alonzo' (2002); 'Jiang, Leroueil, and Konrad' (2004); 'Richefeu, Youssoufi, and 



Radjai, ( ,2006^ ; ,Shamy and Groger, ( ,2008^ ; ,Scholtes, Chareyre, Nicot, and Darve 
( |2009| )). The DEM is essentially a Lagrangian (mesh-free) technique where each 
particle of the material is a sphere identified by its own mass, radius and moment 
of inertia. For every time step of the computation, interaction forces between par- 
ticles, and consequently resulting forces acting on each of them, are deduced from 
sphere positions through the interaction law. Newton's second law is then inte- 
grated through an explicit second-order finite difference scheme to compute new 
spheres' positions. 

We present here a 3D discrete element model developed into the YADE-Open DEM 
platform ( [Kozicki and Donzel ( |2008| )), where effects of capillarity have been imple- 
mented in both terms of force and water retention through the resolution of the 
Young-Laplace equation as presented hereafter. In the proposed model, each in- 
teraction can thus be described as the sum of a dry contact repulsive force and of 
an attractive capillary force resulting from the presence of liquid bridges between 
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Un<0 
(a) ^-_-^ (b) 

Figure 1: Contact interaction: (a) model and (b) overlap. 




(a) 



Un 




Figure 2: Contact friction law: (a) normal force and (b) tangential force. 



grams. 



2.1 Contact friction model 



The elastic-plastic model for contact friction is presented in Fig.[TJa). Although the 
formulation is quite simplistic compared to more realistic Hertz-Mindlin solutions, 



it has been proved to give accurate results ( [DiRenzo and DiMaio| ( |2004| )), with a 
significant computational efficiency. 

Particles are considered to be rigid, but can overlap as shown in Fig. [TJb). This 
overlap accounts for the deformation induced by contacts forces. A linear elastic 
law provides the contact force as a function of the relative displacement between 
two interacting grains (see Fig. |2]). A normal stiffness K^ is defined to relate the 
normal force F^ to the intergranular normal distance Un such as: 



KnUn, if Un^O 



(1) 



"" ^0, if Un>0 
and a tangential stiffness Kt allows to deduce the shear force Ft induced by the 
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(b) 



Figure 3: Illustration of a liquid bridge between two particles of unequal sizes: (a) 
global geometry, (b) details of the bridge. 



incremental tangential relative displacement dUt'. 
dFt = -KtdUt 



(2) 



Kn and Kt are dependent functions of the interacting particle radii R\ and R2 and of 
a characteristic modulus E of the material such as: 



^ _ 9 E.RI.R2 

Kt — a.Kn 



(3) 



with a a fixed parameter. This definition results in a constant ratio between E and 
the effective bulk modulus of the packing, whatever the size of the particles. 



Shear and normal forces are finally related by a slip Coulomb model such that 
pmax ^ —iiFn, whcrc /i is the contact friction coefficient defined as /i = tan{(pc), 
with (pc the intergranular friction angle. 

2.2 Capillary model 

For simplicity, we assume that capillary water inside the sample is solely composed 
of interparticular independent menisci as defined in the pendular state (Fig. [3]). 

Their exact shape between spherical bodies is defined by the Young-Laplace equa- 
tion, which relates the pressure difference Au = Ug — ui across the gas-liquid inter- 
face to the radius of curvature C of the bridge surface and to the surface tension y 
of the liquid phase: 



Aw ■■ 



1 
C 



(4) 
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In the Cartesian coordinates of Fig. [3j C can be formulated as a function of the 
profile y{x) of the liquid-gas interface curve, the x axis coinciding with the axis of 
symmetry of the bridge, passing through the centers of the bonded spheres. Fol- 
lowing the Young-Laplace equation (Eq. |4]), the profile of the liquid bridge is thus 
related to the capillary pressure Aw through the following non-linear differential 
equation: 

^a^,.-(.r^-^'-±m-v'(.^=o (5) 



According to a recent study by [Soulie, Cherblanc, Youssoufi, and Saix| ( |2006| ), 



the corresponding meniscus volume Vm and intergranular distance Un can be ob- 
tained by considering the x-coordinates (x^i and Xc2) of the three-phases contact 
lines defining the solid- liquid-gas interface such as: 



nJ^^fy^{x)dx-^nRl{l-acos{xci))^{2 + acos{xci)) 



(6) 



■^7r/^2(l — cicos{xc2))^ {^2. + acos{xc2)) 
and 

Un—Rii^ — acos{xc2)) +-^c2+^i(l — acos{xci)) —Xci (7) 

The capillary force due to the liquid bridge can then be calculated at the profile 
apex yo according to the 'gorge method' ( |Hotta, Takeda, and Ionya| ( |1974| )) and 
consists of a contribution of both Aw and 7: 

Fcap^27ly^y+7lyl^u (8) 

The relation between the capillary pressure and the configuration of the capillary 
doublet is thus described by a system of non-linear coupled equations (Eq.[5j Eq.|6j 
Eq.[7j Eq.[8]) where local geometry ([/„) and meniscus volume (Vm) arise as a result 
of the solved system. 

In order to account for capillarity in the YADE-Open DEM code, an interpola- 
tion scheme on a set of discrete solutions of the Young-Laplace equation has been 
developed so as to link directly Aw to Fcap and Vm for a given grain-pair configura- 
tion {Ri, R2, Un). This numerical procedure results in a suction-controlled model 
where, at every time-step during the simulation, capillary forces and water volumes 
are computed based upon the local geometry and the imposed capillary pressure. 
Fig. |4] presents the evolution of the capillary force with the relative displacement 
between two interacting grains. 

It is to be noted that capillary forces act exclusively in the axial direction of the 
interaction (the normal to the tangential contact plane), and that Fcap is maximum 
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cap 



Figure 4: Evolution of the capillary force with the relative displacement between 
interacting grains for a given capilary pressure Aw: a meniscus can form for Un < 

U creation ^^^ breaks for Ufi > U rupture • 



for grains in contact {Un < 0). Umpture is the debonding distance corresponding to 
the minimum Un value from which the Young-Laplace equation has no solution. 
An hydraulic hysteresis can be accounted for by defining Ucreation as the distance 
from which a liquid bridge can form between two interacting grains. For instance, 
as a drying scenario results in a desaturation of the material, liquid bridges can exist 
at stretched contacts as long as the debonding distance Un is less than the critical 
separation distance Umpture (Fig-Ra)). On the other hand, a wetting scenario only 
allows appearance of menisci between contacting or close grains (Fig. [5|b)) such 
as it would occur during a capillary condensation when the relative humidity of the 
surrounding air is increased. The choice was made here to allow bridge appearance 
for strict contact {Ucreation = 0), thus neglecting the possible effect of adsorbed 
water. In the same way, once they break, menisci can reform when particles come 
strictly in contact. 

Following these assumptions, wetting and drying scenarios can be taken into con- 
sideration for the initial state of the material, keeping in mind that, in our simula- 
tions, the contact angle 6 (Fig. [3jb)) which defines the wetting of the material is 
set to 0° for simplicity, assuming a perfect wetting of the material and occulting 
therefore its possible effect on local hysteresis mechanisms. 



2.3 Interparticle behaviour 

The interaction laws involve normal repulsion and Coulomb friction at contact as 
well as capillary adhesion as presented respectively in Fig. [2] and Fig.|4j 
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Figure 5: Hydraulic hysteresis history accounted for in the code: assembly of par- 
ticles resulting from (a) a drying scenario, (b) a wetting scenario. 



Capillary forces are considered constant for the range of elastic contact deformation 
(Un < 0), assuming a low deformability of particles during simulations. To pull two 
bonded spheres apart, an increasing external tensile force is needed. This force has 
to balanced the sum of the repulsive contact force F^ and the adhesive capillary 
one Fcap which tends to decrease with the interparticle distance C/^. The tensile 
process is stable as long as Un is less than the debonding distance Umpture which 
will depend on both the doublet configuration {Ri, R2) and the capillary pressure 
Aw. As presented in Fig. 6l the particle size ratio ^ = f^ has a great influence on 
both the adhesive force and the liquid bridge volume. For instance, for a given Au 
value, the greater the ratio r is, the smaller the bridge volume Vm is, resulting in less 
significant interacting distances Umpture before rupture as well as in diminishing 
force intensities. 

Concerning now the influence of the capillary pressure (Fig. |7]), it has to be noted 
that, for a given doublet configuration (r = 1), if higher values of Aw result in 
smaller liquid bridge volumes, capillary force intensity at contact {Un = 0) is quite 
constant whatever Aw. At the particle scale, capillary pressure, through its effect on 
menisci volumes, finally influences more the debonding distance than the resulting 
adhesive force. 

To sum up, in addition to the great influence of particle size distribution on their in- 
tensity level, capillary effects between interacting grains are strongly driven by the 
amount of liquid involved in their behaviour. Such a complexity leads to numerous 
uncertainties when considering granular materials such as soils. The next section 
is therefore an attempt to provide some clarifications on capillarity consequences 
at the scale of a numerical granular assembly. 
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Figure 6: Effect of the particle size ratio r = |^ on (a) the normalized capillary 
force F^^p = ikRy ^^^ ^^^ ^^^ normalized meniscus volume V^ = ^ as a function 
of the separation distance f/^ = ^ for Aw = 1000 kPa. 
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Figure 7: Effect of the capillary pressure Aw on (a) the normalized capillary force 
^cap ^ 27iRy ^^^ ^^^ ^^^ normalized meniscus volume V^ = ^ as a function of the 
separation distance C/^ = ^ for r = |^ = 1. 
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Table 1: Input parameters 



Global Modulus 
E (MPa) 


k, 

kn 

a 


Friction angle 
(Pc (deg.) 


50 


0.5 


30 




Figure 8: The discrete element model. 



3 Macroscopical behaviour 

In order to investigate macroscopic consequences of interparticle capillary menisci, 
numerical simulations were conducted on a polydisperse assembly composed of 
10,000 spheres (Fig.[8]) with a uniform grain size distribution ranging from 0.035 mm 
to 0.07 mm. The specimen is contained inside a box made of six rigid frictionless 
boundary walls for which positions at each time steps are defined from the pre- 
scribed loading program. Tab. [T] provides a summary of the input parameters used 
in the simulations, referring to Eq.[T} Eq.[2]and Eq.[3] 

The sample is prepared in a dry configuration by an isotropic-compaction technique 
which ensures the initial homogeneity of the packing ( [Scholtes, Chareyre, Nicot, 



and Darvel ( |2009| )). The final porosity of the sample depends on the intergranular 
friction angle (pc defined during the compaction stage. The smaller (pc, the denser 
the specimen. Here, (Pc was fixed to 1°, leading to a dense specimen with a porosity 
of about 0.39. (Pc was then set to 30° for all the simulations. Starting from the 
stabilized isotropic configuration, capillary pressure Aw is then switched on, letting 
menisci to be formed between all possible interacting grains following the chosen 
drying or wetting scenario (stretched or contacting menisci respectively). It is to be 
noted that the suction controlled scheme ensures menisci to be homogeneously dis- 
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tributed inside the medium in accordance the thermodynamic equilibrium between 
both liquid and gas phases. 

3. 1 Numerical results 

Several loading paths were applied to the specimen in order to assess both its hy- 
draulic and mechanical properties. 

3.1.1 Hydraulic properties 

Determining hydraulic properties of unsaturated soils is of great interest because 
of the strong hydro-mechanical coupling consequences on the overall behaviour. 
Herein, the proposed model allows to determine the water content of an assembly 
as a result of the summation from all menisci volumes contained inside the medium. 
Indeed, the degree of saturation Sr of the numerical sample is obtained by: 

5r=100*^ = 100* ^•"=^ "" (9) 

^v ^sample ~ ^grains 

with V/ the total liquid volume and Vy the void volume resulting from the grains 
arrangement inside the cubic box. Vy is equal to the difference between the to- 
tal volume of the sample Vsampie minus the volume occupied by all the grains g 
forming the assembly (Vgrains = Lali f^^g)- ^Y varying the capillary pressure Aw 
inside a given sample, one can therefore construct its soil water characteristic curve 
(SWCC). 

However, due to the pendular assumption concerning liquid phase modelling (in- 
dependent interparticular liquid bridges), the SWCC cannot be computed for high 
degrees of saturation Sr. Indeed, below a given capillary pressure, liquid bridges 
become big enough to allow their overlap with neighbours. As this overlap leads 
to menisci fusion associated with a complete fluid reorganisation which makes ob- 
solete the pendular assumption, a numerical procedure has been developed in order 
to identify menisci superpositions through the definition of their geometrical con- 
figuration. As defined in Fig.[3j the wetting of a liquid bridge over the grain surface 
is defined by its filling angles 5i and 82. It is therefore possible to determine if, 
for a given grain g of an assembly, neighbouring menisci meet through their wetted 
surfaces. 

For instance (see Fig. M), if the angle OtOgOj formed between three interacting 
grains g, i and j is smaller than the sum {5gi + 5gj) corresponding to the surface of 
the grain wetted by menisci rrig^i, and rrigj linking g with / and g with j respectively, 
an overlap occurs and fusion can therefore process. Obviously, the occurence of 
menisci fusion strongly depends on the microstructural arrangement and on the 
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(a) 



(b) 



Figure 9: illustration of the numerical procedure for menisci fusion: (a) definition 
of the local geometry, (b) identification of merged menisci: rugj, rug^k and m^ / 
overlap on grain g surface. 



porosity of the medium. The range of saturation that corresponds to the pendular 
regime cannot therefore be infered for all materials and has to be checked for each 
particular case. 

Fig. [To] presents the SWCC obtained for the predefined numerical assembly in the 
range of its pendular regime. For this particular case, our finding is that the inde- 
pendence of menisci is strictly ensured for Sr < 12%, whereas overlappings are 
negligible until 5r = 15 %. 

One can note that, through the possible local mechanisms of formation and break- 
age of menisci, an hydraulic hysteresis arises at the macroscopic level, depending 
on the wetting history of the material. This macroscopical hysteresis is directly 
linked to the liquid bridge density inside the sample. Indeed, due to the possibility 
of stretched liquid bridges between separated grains for a drying scenario, menisci 
are numerous inside the medium and the corresponding degrees of saturation are 
therefore larger for the simulated dried material than for the wetted one. Fig. [TT 
illustrates this point through the evolution of the average number of liquid bridge 
per particle Kyn during both the two scenarios. 

It is remarkable that both curves present two distinct stages: 



- For Aw > 3 MPa, menisci distribution is driven by the particle gradation. 
Menisci progressively disappear with the increase in Aw due to the persis- 
tence of capillary water between the smallest particles. For instance, this is 
the reason why clay needs higher capillary pressures than sand to be fully 
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Figure 10: Soil Water Characteristic Curve generated with the proposed discrete 
element model. 
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Figure 12: Simulation of drained triaxial compressions on the dry specimen under 
5, 10 and 20 kPa confining pressures. 



desaturated. 

- For Aw < 3 MPa, menisci distribution is driven by the capillary pressure ac- 
cording to the wetting history of the material. During the drying scenario, 
the increase of Aw induces smallest debonding distances (Fig. [7]) and there- 
fore less stretched bridges inside the sample. During the wetting scenario, as 
menisci are only present at contacts, they keep constant in number whereas 
their respective volumes increase with the decrease of Aw. 

At the particle scale, the simulated hysteresis can be illustrated through Fig.[5]where 
one can see that the darkest particle has more neighbours (K^ = 5) after a drying 
sequence, than after a wetting one (K^ = 2). This phenomenon, also denoted as the 
ink-bottle effect, is probably the most significant hysteretic mechanisms involved in 
unsaturated soil mechanics. The consequences of this property on the mechanical 
behaviour of granular assemblies is investigated in section 3.2. 

3.1.2 Stress-Strain response under triaxial loadings 

To determine the effect of interparticular capillary menisci on the mechanical prop- 
erties of a granular material, we performed a series of triaxial compression loadings 
(<72 = (73, £1 > 0) on the numerical assembly under 5, 10 and 20 kPa confining pres- 
sures. 



As a reference, the response of the dry specimen is plotted in Fig.[T2j The assembly 
presents a behaviour which is qualitatively close to that of a dense granular material 



with strain softening associated with dilatancy. As shown in Fig. [13} the discrete 
element model behaviour is well described through the Mohr-Coulomb criterion: 



T = G tan^ + c 



(10) 
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Figure 13: Failure envelop for the dry sample. 



with a null apparent cohesion c and an internal friction angle <t> of about 28°. 

In order to cover the entire pendular regime, numerical unsaturated triaxial tests 
were performed at different capillary pressures corresponding to degrees of satura- 
tion ranging from to 10 % along the wetting path. Tab. [2] presents the correspon- 
dence between Aw and initial values of Sr for the considered specimen. 

Table 2: Capillary pressures and corresponding initial degrees of saturation of the 
discrete element model 



Capillary pressure Am [kPa] 


5000 


3000 


80 


50 


30 


20 


Degree of saturation Sr[%] 


0.001 


0.01 


1 


2.5 


5.5 


10 



Fig. [14] displays the responses of the unsaturated samples compared to the dry one 
for triaxial compressions applied under a confining pressure of 10 kPa. For clarity 
reason, only four of the six tests are plotted. As expected, the shear strength is 
greater for unsaturated materials and depends on the degree of saturation. The 
more wetted the sample is, the higher the deviatoric strength is, with a trend for 
dilatancy more pronounced than in the dry case, suggesting a more interlocked 
structure. To sum up all simulated tests, the corresponding failure envelops as well 
as the evolution of the apparent cohesion c with Sr are plotted in Fig. [15] 

As a first observation, one can see that the internal friction angle <t> is independent 
on the degree of saturation, whatever the water content level. This probably re- 
sults from the assumption of null effect of capillarity on contact friction, but it is 
remarkable that such a local property arises at the scale of an assembly. |Richefeu, 
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Figure 14: Simulation of triaxial compressions at different capillary pressure levels 
under a 10 kPa confining pressure. 
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Figure 15: (a) Morh-Coulomb failure envelops and (b) apparent cohesion c as a 
function of the degree of saturation Sr. 
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Figure 16: Mean capillary force < Fcap > as a function of the degree of saturation 
Sr. 



Youssoufi, and Radjai'l ([2006*) came to the same conclusion from laboratory ex- 
periments on glass beads assemblies with pure water, providing therefore a good 
prediction to the idealized model. However, real granular materials such as soil 
would certainly not highlight such an independence between local and global fric- 
tion properties due to the complex combined effects of surface roughness and in 
situ liquid wetting features. 

Second, it is clear that the cohesion varies significantly from dry to unsaturated 
states, c increases non-linearly till a maximum value for the higher degrees of 
saturation that corresponds to the upper limit of the pendular state. This is fairly 
typical based on synthesis of reported experimental data on sand (Lu, Bail in, and 
Tan| ( |200 7)) or glass beads ( [Richefeu, Youssoufi, and Radjai'l ( [2006) )). Nevertheless, 
this is quite remarkable in regards to the interparticle behaviour, since capillary 
force intensity at contact is not such dependent on capillary pressure (Fig. [7]). For 



instance, computing the mean capillary force < Fcap >= 






in the specimen 



for all tested degrees of saturation (Fig. [16]) indicates that capillary forces tends to 
be greater, on average, for higher capillary pressure levels (and therefore for smaller 
saturation degrees). 

This paradox is clearly due to capillary pressure influence on the debonding dis- 
tance of menisci as presented in Fig. [7] Indeed, strong internal rearrangements 
occur during the loading, leading to a redistribution of the liquid inside the mi- 
crostructure which is not obvious at the macroscopic scale. How this redistribution 
occurs and how it influences the shear strength of the material is strongly linked to 
the elaborated properties of grains when they interact in an assembly. A microme- 
chanical investigation is therefore needed in order to identify effects induced by 
structural evolution of the medium on the liquid bridge distribution. 
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Figure 17: (a) Evolution of the average number of menisci per particle K^n during a 
triaxial compression for several initial degree of saturation Sr, (b) K^ values at the 
q-peak as a function of Sr. 



3.2 Microscopic analysis 

Fig. [TTfa) shows the evolution of the average number of menisci per particle K^^ 
for different degrees of saturation during the unsaturated triaxial compressions pre- 



sented before (Fig. [14]). 

Due to particle rearrangements during the loading, liquid bridge distribution inside 
the material changes with deformations, depending on the defined capillary pres- 
sure level. For instance, despite the fact that the initial liquid bridge densities are 
identical for tests corresponding to degrees of saturation of 0.01, 2.5, 10 %, their 
evolutions during compression are totally different due to the local response of liq- 
uid bonds to loading. Typically, as a result of the local behaviour of menisci (see 
Fig. [7]), liquid bridges tend to persist more with the medium dilatancy for low val- 
ues of Am (higher Sr), due to larger debonding distance U rupture- As discussed in 



Scholtes, Chareyre, Nicot, and Darve| ( |2009| ) as well as in |Scholtes, Hicher, Nicot, 
Chareyre, and Darvel ( |200 9'), the microstructure strongly influences liquid distribu- 
tion inside the medium. The initial bridge density number appears therefore not 
appropriate to accurately evaluate the shear strength of the material as it was sug- 
gested by |Richefeu, Youssoufi, and Radjai| ( [2006| ), in the sense that failure does not 
occur for an isotropic configuration of the material, but after some deformations 
leading to significant changes inside the fabric due to the induced anisotropy. Co- 
hesion, as well as the internal friction angle, is dependent on space direction, and 
has to be evaluated at the yield state. Indeed, by plotting Kyn values corresponding 
to the q — peak (Fig. [T7|b)), one can recover the shear strength hierarchy induced 
by capillarity depending on the degree of saturation. 
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Figure 18: (a) stress-strain relationships and (b) evolution of K^ during a triaxial 
compression on a sample subjected to a drying and a wetting scenario respectively 
under the same capillary pressure Au=20 kPa. 



It is important to note that, in the end, whatever the capillary pressure level or the 
degree of saturation, the pertinent parameter which determines the shear strength of 
a partially saturated granular material is the liquid bridge density inside the medium 
at the yield state. Obviously, the evolution of this density with deformations is 
totally driven by the capillary pressure level, and, thereafter, by the water content. 

Another way to illustrate the primacy of the liquid bridge density on shear strength 
properties of a wet granular material is to compare its behaviour for both a drying 
and a wetting scenario. As suggested in section 3.1.1, a given sample can contain 
different numbers of capillary bonds for a defined Aw value (Fig. [TT]) due to the 
possible hydraulic hysteresis. 

Fig. [TSja) presents the stress-strain relationships obtained for a numerical sample 
subjected to a triaxial loading under the same capillary pressure (Aw = 20 kPa), but 
with two different initial configurations resulting respectively from a drying and a 
wetting scenario. Due to numerous liquid bridges Fig.[T8jb), the dried sample has 
a greater shear strength (qpeak=^^ kpa) than the wetted one {qpeak-'^'^ kpa), which 
results in respective cohesions of 5 and 6.5 kPa. It is remarkable that the ratio 
^^wetting = 1.3 is Strongly correlated to the ratio between the corresponding values 

j^Drying rj . 

of Km at the q — peak: wettmg ^ 53 ^ 1-34, confirming here the possibility to link 

f^^-pea ^^ ^^^ overall cohesion of the material. 

Another remarkable feature of the coupled hydro-mechanical process is that in- 
ternal rearrangements caused by the loading tend to conceal the initial difference 
between the two specimens, leading to a common residual state for large deforma- 
tions with the same number of liquid bridges on average. This confirms the strong 
influence of the microstructure on the liquid phase distribution. For instance, if the 
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presence of fluids inside a polyphasic granular material implies some consequences 
on its overall behaviour, the coupling is also active in the other way. Microstruc- 
tural rearrangements of the solid skeleton lead to strong modifications of the fluid 
phases and both aspects induce complex mechanisms at the macroscopic scale. 

4 Conclusion 

A micromechanical computational model for the analysis of wet granular soils in 
the pendular regime has been proposed. Capillary mechanisms are described at the 
contact scale based upon the capillary theory in both terms of interparticle adhe- 
sive force and water retention. Capillary menisci are distributed inside the medium 
according to the thermodynamic equilibrium between liquid and gas phases. More- 
over, an hydraulic hysteresis is accounted for based on the possible mechanisms of 
formation and breakage of liquid bridges during wetting and drying phases. 

Triaxial compression test simulations were performed on a granular assembly un- 
der several confining pressures for dry and partially saturated conditions in order to 
analyze effects of local capillary menisci on macroscale cohesion and friction prop- 
erties. The results confirm that capillary-induced attractive forces and hydraulic 
hysteresis play a significant role in the shear strength of granular materials. By 
increasing normal forces at contact, capillary menisci contribute to the apparent 
cohesion of the material and enhances its stiffness. A remarkable aspect is that the 
shear strength of an assembly is basically controlled by the liquid bridge density 
and not by the capillary pressure level. A strong correlation is highlighted between 
the average number of liquid bridges at the yield state and the induced cohesion. 
On the other hand, the liquid bridge distribution inside the medium is controlled by 
the capillary pressure which leads to counter-intuitive response at the scale of an 
assembly. In a separate way to capillary forces. Coulomb cohesion increases with 
water content and tends to a maximum value for the upper limit of the pendular 
regime following the liquid bridge density. In additon, the results put in evidence 
the effect of the skeleton induced anisotropy on liquid phase properties which is 
generally occulted when considering hydro-mechanical modelling. Solid fabric 
evolution clearly has a strong influence on the liquid distribution. 

The proposed model is able to simulate the macroscopic response of wet granu- 
lar materials and revealed a number of outstanding micromechanical mechanisms 
and response patterns consistent with experimental data. The multi-scale approach 
presented here appears to be a pertinent complementary tool for the study of unsat- 
urated soil mechanics. A fundamental extension of this work would be to extend 
the interaction model to higher degrees of saturation with, for example, a validation 
through tomographic imaging. 
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